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Angle resolved photoemission experiments by Kanigel, et al [Nature Physics 2, 447 (2006)] have 
made a remarkable observation that low energy electronic excitations in the normal state of un- 
derdoped cuprate superconductors exist on open "Fermi arcs" instead of a closed Fermi surface. 
These arcs shrink upon cooling, with the arc length appearing to extrapolate to nodal points at zero 
temperature. We show that this striking non-Fermi liquid behavior could result from the under- 
doped normal state above Tc lying in the vicinity of a quantum phase transition between a d-wave 
superconductor and a correlated insulating phase. 
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Introduction. — The normal state of the underdoped 
high temperature superconductors (SCs) has long been 
argued to lie outside the paradigm of Fermi liquid the- 
ory due to strong electron correlations. The most re- 
markable manifestation of this appears in recent angle re- 
solved photoemission spectroscopy (ARPES) studies [l| 
above the superconducting transition temperature (Tc) 
in imderdoped (UD) Bi2Sr2CaCu208+5. These experi- 
ments shovif that the spectral vifcight of low energy single 
particle excitations is not present on a closed Fermi sur- 
face (FS) in momentum space. Instead, it forms open 
"Fermi arcs" [ij of length - T/T*, where T* is the pseu- 
dogap temperature, indicating a complete breakdown of 
the Fermi surface between Tc and T*. One can envision 
two very different causes for such NFL behavior. 

(i) The ARPES experiments could be probing the finite 
temperature properties of some exotic non-SC ground 
state which competes with d-wave superconductivity at 
low doping. Since the Fermi arcs appear to extrapolate 
to nodal points enclosing zero volume [l| as T — > 0, this 
"underlying normal state" must be a non-Fermi liquid 
as it violates Luttinger's theorem [3, Varma and 
Zhu have postulated a time-reversal symmetry break- 
ing phase with circulating currents and argued that cur- 
rent fluctuations in this phase can lead to Fermi arcs. 
Other proposals for candidate underlying normal states 
include the d-density wave Q and staggered flux phases 
Q- With varying doping, these two states support hole 
pocket Fermi surfaces with an area proportional to the 
doping, although the spectral weight around this pocket 
is anisotropic and has been argued to lead to an arc-like 
feature. Another concern with these two states is that 
the "arcs" in such models do not appear to extrapolate 
into nodal points as T — > 0. One more recent proposal 
for arcs involves inducing an additional nematic liquid 
crystalline order in the SC state at low doping Q but so 
far it has not addressed the normal state above Tc- 

(ii) A second possibility is to examine whether order 
parameter fluctuations of the d-wave SC or fluctuations 
of the short range antiferromagnetic order above Tc could 
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FIG. 1: (color online) Temperature evolution of Fermi arcs, 
regions of significant j4(k,tj — 0) in the first quadrant of the 
Brillouin zone. 

lead to Fermi arc physics. This has been examined ear- 
lier using various approaches although these works 
have not focussed on the temperature scaling of the arc 
length. A distinct way to study thermal order parame- 
ter fluctuations is to consider the effect of a dilute gas 
of classical vortices [l^l above Tc- Alternatively, at low 
dopings and temperatures, quantum fluctuations of the 
order parameter are also relevant; such quantum fluctu- 
ations could be particularly important if one is close to 
a quantum phase transition out of the superconducting 
state into a correlated insulating state. 

In this paper we will explore this last route. Namely, 
we study the possibility that the Fermi arcs in the under- 
doped normal state arise from proximity to a quantum 
phase transition between a d-wave SC and a correlated 
insulating state at low doping. Fig. 1 illustrates one of 
the key results of this paper — the computed zero en- 
ergy spectral function for the model which we study in 
this paper shows temperature dependent Fermi arcs in 
the normal state as seen in ARPES. 

There are several experimental motivations to explore 
this route which we follow. First, the observed gap struc- 
ture in the normal state appears to extrapolate smoothly 
with lowering temperature to a SC gap with a pure d- 
wave form over a range of dopings. This suggests that any 
"underlying normal state" scenario has to be speciflcally 
chosen to extrapolate into the excitation spectrum of the 
d-wave state over a whole range of dopings. This seems, 
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to us, to be slightly unnatural although we cannot rule it 
out. It appears more natural to assume that the normal 
state has fluctuating d-wave correlations. Second, the arc 
length seen in the experiments appears to scale linearly 
with temperature and independent of most microscopic 
details. Such a simple T-dependence suggests proximity 
to a quantum phase transition near which the temper- 
ature naturally arises as the only energy scale. Finally, 
NFL behavior has been observed in a variety of other 
materials [ill [T^ near field or pressure tuned quantum 
phase transitions. Given this ubiquity of NFL behavior in 
the vicinity of quantum phase transitions, it is worth ex- 
amining this scenario also for the underdopcd cupratcs. 
Our proposal receives support from penetration depth 
measurements [l3 on YBa2Cu306+2; which appear to be 
probing the low temperature manifestations of a 3D SC- 
insulator quantum phase transition [H, 15|. Recent tun- 
neling experiments [lit other cuprates also show that 
the SC state at low hole density is destroyed by charge 
localization leading to such a SC-insulator transition. 

SC-insulator transition, critical charge fluctuations 
and the electron spectral function. — The underdoped 
cuprate SCs are doped correlated insulators, and the 
strongly correlated Hubbard model is thought to be a 



possible minimal description of this system 17|. At 
large U, this leads to an effective tJ model with a ki- 
netic energy for the doped holes and a superexchange 
interaction between the spins. In order to calculate the 
electron spectral function in such a system, we appeal 
to slave-boson mean field theory [B, S, [11, ^ which 
decomposes the electron operator into a product of a 
spin-1/2 neutral spinon and a spinless charge-1 boson, 
cj.(r, t) ~ /|(r, r)^'(r, r). These two particles are cou- 
pled by gauge fields enforcing the non-double occupancy 
constraint. At mean field level [l7, 18|, the spinous and 
charges are decoupled since gauge fluctuations are ig- 
nored — the SC state then corresponds to d-wave pairing 
of spinous and condensation of the bosons. Despite this 
drastic approximation, the SC correlations obtained in 
this theory are in remarkable qualitative agreement with 
a careful treatment using Gutzwiller projected wavefunc- 
tions. This wavefunction approach works directly with 
electronic degrees of freedom, and has been shown to 
quantitatively explain a variety of photocmission and op- 
tical conductivity data 2l| in the SC state. At finite tem- 
perature, gauge fluctuations around the mean field state 
are more important; however their effect has been shown 
to be mitigated in the presence of critical boson fluctu- 
ations [i^l (which we are interested in here), although 
the gauge fluctuation corrections to the results presented 
here remains an important problem for future work. For 
now, we will therefore view slave boson mean field theory 
as a crude but qualitatively reasonable way to proceed, 
in order to see what aspects of the phenomenology of the 
underdopcd cupratcs can be possibly captured. 

With undcrdoping, the SC undergoes a quantum phase 



transition into a correlated insulating state which, as 
suggested by tunneling experiments [l^, has glassy 
bond-centered charge order which could arise from the 
Coulomb interactions between the bosons. However, in 
order to expose the physics of the Fermi arcs we choose 
to study a simpler model of a SC-insulator transition, 
where the boson field ^'(r, r) is described by a Ginzburg 
Landau action 



S^^j drd^v (|a,*p + c2|V*p+m2|*|2 + ||5-| 
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Note the SC is three-dimensional and anisotropic, 
c^jV^p is the short hand notation for cIlVn^Pp + 

c^l^z'I'P , with Cz <^ C||. For > 0, the system is in a 
uniform insulating phase with a charge gap. For < 0, 
the SC state with a nonzero (4') is more stable. The mean 
field SC-insulator critical point is at = 0. Model ([T|) is 
at its upper critical dimension and we can drop the quar- 
tic interaction at T = 0. The critical boson Green func- 
tion is then x{Q,^l'n) = (*q,,^„ *Q,<^J = + ^q), 
where Q = (q, Qz ) is the 3D momentum with q being its 
ab-plane component, Qq = c|q^ -|- c^g^, and i/n ~ 2mrT. 
Sitting at m? = and warming up from zero tempera- 
ture, the quartic interaction induces [2^ a T-depcndcnt 
mass R leading to the charge fluctuation propagator 
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(2) 



Leading order perturbation theory in g at 1-loop 
(Hartrec) level yields R — 4? f2Q"'^n(f2Q), where 
n{x) = l/(e^/^ — 1) is the Bose function. At finite tem- 
perature, the interaction also induces a damping term in 
the boson Green function. This thermal relaxation rate 
is however higher order in g in a perturbative treatment, 
and we ignore it here. A description of the thermal re- 
laxation rate over a range of energies from small to large 
values of v/T is nontrivial [l^l; we have so far not pur- 
sued this calculation. 

We assume the spinous are paired in the d-wave 
channel throughout the transition, and described by 
a d-wave Bardeen-Cooper-Schrieffer type Hamiltonian 
with a kinetic energy = —2t{coskx + cosfcj,) — 
At' cos kx cos ky — /i, and a d-wavc pairing gap Ak = 
A (cos kx — cos ky)/2, with the resulting quasiparticle dis- 
persion i?k = ■\/^k + ^k- Support for this approximation 
comes from thermal transport measurements [31 which 
indicate that the thermal conductivity, which is domi- 
nated by nodal d-wavc quasiparticles in the SC state, re- 
mains unchanged when we enter the non- superconducting 
state below a critical doping. Using this, the electron 
Green function can be obtained by convolving its spinon 
and boson parts. 



G(k,ifc„) = ^ 
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where N is the number of sites, 

1 + n(a;Q) - /(i;k-q) , ^(c^q) + /(i;k-q) 



Fi 
F2 



ikn ~ i?k-q - ti^Q ikn ~ -^k-q + Wq ' 

1 + n(u;Q) - /(i;k-q) , n(coQ) + /(iJk-q) 



lfc„ + i?k-q + l^Q 



Here cjq = y'fi^^ + R, f{x) = l/(e^/^ + 1) is the Fermi 

function, and = (1 + ^k/£^k)/2 and v^. ^ 1 — are 
the usual BCS coherence factors. Our work is quite close 
in spirit to the zero temperature calculations of Ref. [2^ . 

After analytic continuation, ikn — > w + zO+ , the imag- 
inary part of Eq. ^ gives the electron spectral function 
A(]<.,uj). The result is easy to understand qualitatively. 
An electron with energy u> and in-plane momentum k is 
made by combining a spinon and a boson. However, the 
spinon Green function has weight at energies ibE'k-q and 
the boson excitations also has weight at ±f2q (as seen 
from the boson Green function in Eq. [2]) which leads to 
four contributions to the electron spectral function. The 
thermal and coherence factors appearing in the spectral 
function appropriately weight the probability of such pro- 
cesses. 

Fermi arcs in A(k,0). — We first argue that Eq. ([3]) 
leads to Fermi arcs centered around the nodes of the d- 
wave SC dispersion if we focus on the zero energy spectral 
function A(k, 0). In order to make the qualitative argu- 
ment, let us work at very low temperatures such that 
T -C A,C||A, c_lA. In this low temperature regime, the 
interaction induced mass R ^ T^. Thus, the charac- 
teristic boson fluctuations will have very small in-plane 
momenta |q| ~ (T/c||) ^ A, small overall momentum 
Ql ^ A, and energies ~ T. In order to obtain a zero 
energy electron with a certain momentum k, we have to 
use the energy of this thermal boson fluctuation to ex- 
cite a spinon quasiparticle, namely spinon quasiparticles 
with binding energy ~ T. The zero energy electron will 
have nearly the same in-plane momentum as the spinon 
quasiparticle (since the boson carries very little momen- 
tum). Given these constraints, and the fact that spinon 
quasiparticles with energy ~ T lie on banana shaped arcs 
of size ~ T/A around the d-wave nodes of the SC, it is 
clear that the zero energy electrons will also lie on "Fermi 
arcs" of size T/A around the nodes of the underlying d- 
wave SC. Since T* ^ A in the cuprates, this is consistent 
with (T/T*) scaling of Fermi arcs reported in Ref. [l[. 

The above rough argument relied on working at very 
low temperature. We next proceed to a numerical eval- 
uation of the electron spectral function A(k, 0) in order 
to justify this argument and also to study the evolution 
upto temperatures T ^ A. We have numerically evalu- 
ated A(k, 0) using Eq.Q for temperatures of interest to 
the ARPES experiments (T = O.lA - 0.9A). We work 
at doping S = 0.05 and choose t' /t = -0.25, A/t = 0.15 
for the spinous. Also we set C||A = 3t :s> A, Cz = C||/300, 
and use a Lorentzian broadening of O.Olt in Eq.®. The 
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FIG. 2: (color online) Estimate for the arc length A<j) from 
the zero energy spectral weight A{<f), 0) for angles (j) along the 
FS (defined in Fig. 1). (a) Arc length defined as angular 
distance between points of maximum A{<p, 0) along the Fermi 
surface, (b) Arc length defined as the region with half the 
total spectral weight along the FS. See text for details. 

results for g = 0.1 are plotted in Fig. 1 for < kx,ky < tt, 
qualitatively similar results are found for other values of 
g. We see from Fig.l(a-c) that A(k, 0) forms open arcs 
in momentum space with negligible (although nonzero) 
weight at other momenta on the Fermi surface. 

Size of the Fermi arcs. — In order to measure the 
size of the Fermi arcs, we analyze the zero energy mo- 
mentum distribution curve along the spinon Fermi sur- 
face, A{(p,0), where is the angle from the (0,0)-(7r, tt) 
diagonal as shown in Fig. 1(a). Going from the nodal 
towards the antinodal direction along the spinon Fermi 
surface, A(0, 0) first reaches its peak value at (j)max and 
then drops to a minimum (background) value of Abg at 
(j> ~ 7r/4. We use Acj) = 2(f)max, i-c. the angular span 
between two spectral peaks along the Fermi surface, as 
a first measure of the arc size. Its approximately linear 
T-dependence is shown in Fig. 2(a). 

An alternative measure for the arc size is given by 
A(j) = 2(ph, where the region < (j) < accounts for 
half the total integrated spectral weight, J^'^ d(f>[A{(f), 0) — 

Abg] = i /J"^" d(j)[A{(j), 0)~Abg]. Acj) defined in this way IS 
shown in Fig. 2(b). We see that both definitions yield the 
similar qualitative result that A(j) scales linearly with T. 
This is consistent with our earlier arguments. (The quan- 
titative agreement between the two definitions above is 
sensitive to the details of what percentage of the spectral 
weight is chosen in the second definition.) 

While the above two definitions of the arc length are 
accessible theoretically, they are harder to measure in ex- 
periments due, atleast, to momentum dependent matrix 
elements. We therefore turn to results for the EDC at 
fixed momentum. 

Energy distribution curve (EDC) lineshape. — Fig. 3 
shows the EDC, namely A{4),uj), plotted over a range of 
temperatures and angles along the Fermi surface. At all 
angles, we find that there is a filling-in of low energy spec- 
tral weight as the temperature is raised. In our model, 
this spectral weight growth is facilitated by the tempera- 
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FIG. 3: (color online) Energy distribution curve A{(j},LL!) for 
various angles (^ along the FS, and for increasing tempera- 
ture showing the filling in of the low energy spectral weight. 
For angles closer to the node as in (a), the filling in of the 
spectral weight occurs at lower temperatures consistent with 
experiments. See text for discussion. 

ture dependent mass R which transfers weight from high 
to low energies. The recovery of zero energy spectral 
weight occurs at lower temperature for near-nodal points 
(smaller cf>) than at angles closer to the antinodal direc- 
tion. This growth of the low energy spectral weight is 
such that the dip-like EDC at low energy evolves into a 
peak-like structure which at sufficiently high temperature 
becomes a dominant feature. 

The above results are broadly consistent with the ex- 
perimental observations While the detailed EDC line- 
shape is not quite like the experimental lineshape, the 
comparison is expected to improve if we take into account 
the following points which have been ignored in our anal- 
ysis, (i) The shift of the higher energy edge-like feature 
to larger uj could be somewhat offset by a temperature- 
dependent spinon pairing gap as the local singlet corre- 
lations weaken with increasing T — this needs a more 
careful self-consistent calculation of the spinon gap. (ii) 
The low energy lineshape may improve if relaxational 
dynamics for the bosons [2^ is taken into account, (iii) 
Our model Eq.Q also ceases to be a good approxima- 
tion at large energies where other inelastic channels, in- 
cluding spin fluctuations, become important, (iv) The 
nonzero energy resolution of the ARPES experiments 
(^ 15 — 30meV) has not been taken into account here. 
Given these limitations, we have not attempted a more 
detailed analysis to extract the arc length from the EDC 
as in the experiments. 

Discussion. — We have explicitly shown, within a phe- 
nomcnological slave-boson approach, that critical charge 
fluctuations near a SC-insulator transition as described 
by Eq. ll]), together with spinon pairing in the d-wave 
channel, gives rise to Fermi arcs with size ~ T/T*. We 
finally present some arguments for how the normal state 
of weakly doped cuprates probed in ARPES experiments 
[H could lie in the quantum critical regime of such a tran- 
sition. Recent tunneling experiments [l6| in the weakly 
doped cuprates indicate that the SC state at low doping 
undergoes a charge localization into an insulating state. 
Such a charge localization transition, where the doped 




FIG. 4: (color online) Tunneling and thermal conductivity 
measurements suggest that the cuprate SCs evolve, with dop- 
ing, from a glassy charge insulator with nodal spinon excita- 
tions into a 3D d-wave SC. In the weakly doped regime above 
Tc, we expect quantum charge fluctuations and Fermi arcs 
similar to that in the simple model studied here. 

holes become immobile, is consistent with a SC-insulator 
between a d-wave SC and an insulating state which re- 
tains d-wave paired spinous. Some hints for the existence 
of d-wave nodal spinous at low temperature comes from 
thermal conductivity measurements in samples with 
Tc = 0, which show a residual k/T very close to that of 
the SC state at higher doping. This scenario for the SC- 
insulator transition is sketched in Fig. 4. Warming up 
above Tc in the low doping regime could thus land us in 
the quantum critical regime of a SC-insulator transition. 
Although we have not specifically studied this transition 
from the SC to a glassy charge-localized insulator with 
d-wave spinous, we believe our results strongly suggest 
that this mechanism of critical charge fluctuations near 
a SC-insulator transition could lead to T-lincar Fermi 
arcs in the underdoped normal state above Tc. It would 
be interesting to look for more direct experimental sig- 
natures of such quantum fluctuations in the underdoped 
normal state. 
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